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Abstract 

We present a new Monte Carlo tool that computes full tree-level matrix elements in 
high-energy physics. The program accepts user-defined models and has no restrictions on 
the process multiplicity. To achieve acceptable performance, CAMORRA evaluates the matrix 
elements in a recursive way by combining off-shell currents. Furthermore, CAMORRA can 
be used to compute amplitudes involving continuous color and helicity final states. 

1 Introduction 

In the present LHC- and Tevatron-physics dominated era, the need for understanding of pro- 
cesses with several hard final-state partons has emerged. As the number of Feynman diagrams 
grows factorially with the final state multiplicity, new recursive methods have been developed 
[H [2] and implemented O U [5]. The recursive method essentially factorizes the amplitude 
into contributing off-shell currents, which are stored in the computer's memory, so that the 
number of operations grows only exponentially. Moreover the method allows a fully numer- 
ical approach, where an initialization run determines which currents can be combined to a 
higher-level one, and the matrix element is evaluated by applying the recursive relations in 
the vertex tree. The input from the physics model is coded in a set of recursive relations, 
which are trivially constructed from the Feynman rules. Hence there is no reason why the 
recursive method should be limited to standard-model amplitudes, especially since many new 
physics signatures are associated with multi-jet events or long decay chains at the LHC. Full 
new-physics matrix elements, treating signal and background on equal footing, seem indis- 
pensable to the BSM phenomenologist's toolbox. 

The purpose of CAMORRJij^ is to provide a stand-alone, user-friendly, modular and model- 
independent implementation of the recursive algorithm. As such, the library can serve as 
a matrix-element computing engine for a Monte-Carlo generator, or be used to reweigh an 

^Abbreviation of CAravaglios-MORetti Recursive Algorithm, freely available at 
http : / / www. nikhef . nl / ~ vdoord 
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Figure 1: Outline of Camorra's basic flow structure. 



existing sample of events. For the latter purpose a Les-Houches event ^ interface is pro- 
vided, which reads events sequentially and automatically builds the subprocess vertex trees. 
To achieve a performance competitive with similar Fortran codes, the library relies heavily 
on template techniques, allowing the compiler to inline the recursive relations and optimize 
the resulting machine code. Furthermore, multiplications by constant, sparse matrices such 
as the Dirac matrices are coded by hand and dynamic binding is avoided wherever possible. 

Together with the core routines, various ready-to-use models are included in the package, 
as well as a flat phase space generator, and generators for color and helicity configurations. 



2 Program Outline 

A rough sketch of the library is given by the flow chart in Fig. [T} All classes and functions 
of the library reside in the namespace Camorra, and all header files will be installed in a 
subdirectory named Camorra, so a program using the library will typically begin as follows: 

#include<Camorra/ CM.algorithm . h> 
using namespace Camorra; 



2.1 Algorithm class 

The main class which distributes a computation is called CM_algorithm. It takes the model 
type and incoming and outgoing particles as template parameters, which are therefore compile- 
time constants and common to all subprocesses. Its constructor and process-addition method 
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take character string arguments, e.g. 

CM_algorithiii<QCD,2,3>myalgo("g,g > t,tbar,g"); 
myalgo.add_process("u,ubar > t,tbar,g"); 

At construction, the algorithm class will first create a (unique) instance of the model, and 
a list of all possible currents for the process multiplicity. This list is static, which means 
that all subprocesses and algorithm instances of the same multiplicity in the program will 
make use of it. Furthermore it decomposes the subprocess character strings and matches 
these with particle names, creating a list of process objects. This list is ordered for efficient 
lookup and double-counted processes are erased using the algorithm's loadO method. Calling 

myalgo . construct_trees () ; 

sets up vertex trees for the subprocesses. This routine consecutively performs the follow- 
ing actions: 

• Choose a final-current particle (by default the first one, but this can be manually con- 
trolled) and set up the zero-level currents. 

• Recursively recombine currents into higher-level ones according the the model's vertex 
content. 

• When the highest-level currents are constructed, attempt to match with the final current. 
Recursively run down the tree, deleting all vertices combining to non-matching highest- 
level currents. 

• Determine Fermi signs and fermion flow reversal flags. 

The algorithm also contains an iterator running over the process list. It is controlled by the 
algorithm's set_process member function, taking a character string argument of the same 
form as before. For faster lookup one may use set_process_f lavours or set_process_pdg_- 
ids, which take integer vectors as arguments. Note that if the process is not found, CAMORRA 
will attempt to insert it into the list and construct a vertex tree on the fly. Once a subpro- 
cess is selected, the user attains access to the phase space instances of the external particles 
by calling myalgo. get_phase_space(i). This interface (see appendix [A|) depends on some 
compile-time data stored in the model class, which will be explained in greater detail in the 
upcoming section. The evaluation of the amplitude proceeds by calling 

myalgo . evaluate () ; 

The evaluation first cleans the list of currents, setting all participating sub-amplitudes to 
zero. Then it runs over its list of vertices, applying the series of recursive relations to the 
off-shell currents. Finally it returns the contraction between the highest-level current and 
final one. The algorithm class also provides a spin-summation version of the evaluation func- 
tion called evaluate_spin_sum() . A helicity summation flag for external particle i is set by 
myalgo . sum_spin(i) , the collective flag is set by calling the sum_spins() member. 
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2.2 LH interface class 



Interfacing with other Monte Carlo tools via the Les Houches standard format is achieved by 
constructing a object of the LH_interface type. This is a daughter class of CM_algorithm, 
with extra methods for reading and writing Les Houches Accord event files. Les Houches in- 
put starts with the initialization block readout via the member function rea.difilename.lhe') , 
followed by the sequential event block readout functions read_event () , f ind.process () and 
input.event () . The first copies the data into a memory buffer and returns a positive boolean 
if the event is compatible with the algorithm, the second searches the process list and occasion- 
ally creates a new tree and the third copies the momenta and colors to the tree. If some masses 
in the event file do not match those in the model file (effective quark masses, for example), 
the member function rescale_phase_space () rescales the rest-frame 3-momenta such that 
the invariant masses match the ones specified in the model. Output to LHE files is trivially 
obtained by calling the write (/i/ename.Z/te) method at initialization and the write_event() 
function subsequently. 



3 Model implementation 

The specialty of the CAMORRA library is the ability to deal with a broad range of user-defined 
models. The current version supports theories containing scalar and vector particles and 
fermions. For the latter both Dirac and Majorana representations are allowed. On the vertex 
side, CAMORRA supports a range of built-in Lorentz and color structures which can be arbitrarily 
combined. The modular structure of the package makes it possible to construct new structures 
without modifying any of the existing code. The vertex multiplicity is however limited to four. 
Effective theories containing higher rank vertices can be included by introducing auxiliary 
fields, which are also supported by the library. Furthermore it should be noted that all 
parameters in a model are passed by reference and hence can be modified at runtime without 
having to invoke any function updating the recursive relations. 



3.1 Compile-time data 

To implement a new physics model, one derives a class from the modeK derived model> class 
template. The relevant data to CAMORRA contained in this class is separated into two groups: 
global, compile-time statements and a constructor consisting of a series of particle and vertex 
additions. The model's header file would for example look as follows, 



class myBSM: public model<myBSM> 
{ 

typedef double value.type; 
typedef Minkowski_type spacetime_type ; 
typedef Pauli.basis Dirac_algebra_type ; 
typedef helicity.type spin_vector_type ; 
typedef colour _f low colour .treatment ; 
static const unsigned diniension=4; 
static const int beam_direction=3; 
static const bool coloured=l ; 
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static const bool continuous_helicities=l ; 
static const bool continuous_colours=0; 
myBSMO ; 



The first line represents the numerical type used throughout the computation of the ma- 
trix elements. One can choose any floating-point type for which the standard library's math 
functions are supported. The second and fifth line determine the space-time metric used 
throughout the program. Choosing Euclidean.type compiles all vector contractions to Eu- 
clidean inner products. However, in the current version the spinor and vector wave functions 
are only defined in four-dimensional Minkowski space. The second line chooses the basis of 
gamma matrices used in CAMORRA. This option is mainly implemented for checking purposes, 
although changing it to Weyl_basis may yield a slight speedup. Together with the beam direc- 
tion definition it gives the algorithm enough information to construct polarization vectors and 
massless spinors. Massive spinor wave functions however require a definition of the spin vector 
type-given here on line three-which represents the choice for the construction of a spin vector 
appearing in the massive spinor formulas (see section 4.1). The boolean coloured should be 
defined and set positive whenever the model contains multiplets of (unbroken) symmetries. If 
it is set false, CAMORRA runs in colorless mode, and no definition of colour_treatment is nec- 
essary. If true, the program needs this information for the treatment of objects transforming 
under the SU{N) adjoint representation. Choosing adjoint treats gluons as octets and all 
color structures accordingly, while for the choice above the gluons are represented as qq pairs 
in color space and all color Feynman rules are expanded in this basis ff]. Finally there are 
the continuous color and helicity fiags, controlling the phase space integration method. These 
values determine the interface to the phase space instances of external particles in a process 
(for more details, see appendix [A]) . 



3.2 Particle insertion 

The final line in the code snippet above declares the constructor of the model. Its implemen- 
tation should consist of a series of statements that define the particle and vertex content of 
the theory. Particle insertion methods are 



a.dd_sca.laLZ<representation>(argument list) ; 
a.dd_f ermioiL<representation> (argument list) ; 
a.dd_vectoz<representation,gauge> (argument list) ; 



Here representation is an optional template parameter type denoting the color degrees of 
freedom of the field. Currently fundainental_rep< SU<N> > or adjoint_rep< SU<N> > are 
supported. Representation types may be chained using the composeo class template to de- 
fine particles with multiple colors. The gauge parameter in the vector inserter function can 
be either Feynman_gauge, unitary _gauge or R_vector_gauge. The argument list consist of 
a character string denoting the name of the particle, pointers to the mass and width and an 
integer denoting the particle data group code. Except for the name, any of these fields may 
be omitted to yield zero, e.g.: 

• add_vector<Feynman_gauge>( "gamma" ,22): insert a massless vector in the Feynman 
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gauge with pdg id 22, 

• add_f ermioiK adjoint_rep< SU<3> > >("'^g" ,&m,&w): insert a (Majorana) fermion 
octet with mass m and width w. 

Particle-anti-particle pairs are inserted with the corresponding functions add.scalars, add_- 
f ermions, add_vectors, which take two names in their argument list, e.g. 

add_ferinions< fundament al_rep< SU<3> > >("b" , "bbar" ,&inb,5) ; 
3.3 Vertex insertion 

After the series of particle definitions, the vertex insertion functions can be used to define 
interactions. These are of the form 

add_vertex< co/or structure, spacetime structure> (particles, couplings) ; 

where the argument list contains three or four particle names and the couplings are pointers 
to complex floating-point numbers. Usually only one coupling is allowed, except for the V-A 
type Lorentz structures which require two (see table [T]). Some enlightening examples are 

• add_vertex<vff>("gamnia" , "e+" , "e-" ,&e): inserts a 7^-type coupling in QED with 
coupling constant stored in the variable e. 

• add_vertex<colour_tensor : : d<inp/ei , 1 ,2> , vf f VA>("Z" ,ubar" , "u" ,&cV, &cA):inserts 
a coupling of the form 

where the values cy and ca have been assigned to the variables cV and cA. In the color 
structure we have substituted f undamental_rep< SU<3> > by triplet. Note that all color 
structures in table [2] reside in the namespace colour .tensor and depend on template 
parameters themselves. For the Kronecker delta d, the first one is the representation 
type of the contracted indices, the second and third denote the fields carrying these 
color indices. 

• The squark-quark-gluino vertex in supersymmetric QCD, 
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has a nontrivial ordering of the color indices, as the default ordering of the SU{N) gen- 
erator colour -tensor : -.Ktriplet ,0 ,1 ,2> assumes the first particle to be in the adjoint 
representation and the second and third in the (anti)-fundamental representations (see 
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Table 1: Available Feynman rules in the package. The third column denotes the types of the particle 
arguments in the ad.d_vertex<?;erte2; class>( . . . ) function, where S denotes a scalar, V a vector and 
F a fermion type. The number of coupling constants is listed in the fourth column, represented by 
Cl, C2, . . . in the last column. 



table [2]). The insertion code for the above vertex looks like (omitting the color tensor 
namespace for brevity) 

add_vertex<T<irzp/ei ,2,0, l>,sffR>("~u_L+" , "ubar" , "~g") ; 

During construction of the vertices, CAMORRA checks if the indices of the particles in the 
argument list match those of the Feynman rule parameter. If this fails, the vertex is omitted 
from the internal list and an error message is issued in the log file. Note however that the 
program checks only if index ranges are in correspondence with the recursive equations, so 
that for example a vector particle appearing in a spinor entry of a vertex inserter will not be 
noticed (at least, in four dimensions). Again there is the option to compose color structures 
for vertices with multiply-colored legs with the template type compose<. . .>. 
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Table 2: Available color structures in the package. The first template argument of the color tensor 
class is either a group type (G) or a representation type (R). The other arguments are distinct integers 
i, j, . . . which should be smaller than the rank of the vertex, and which determine the ordering of 
the multiplet argument list in the third column. In the most right column the adjoint-representation 
indices are labeled a, 6, . . . and other representation indices by a, (5, . . .. 



3.4 Additional features 

Finally there are some extra features in the model base class. First, there is a function that 
defines a particle family, which can serve as an argument for the algorithm class to construct 
a list of subprocesses in one go. For example, the particle content of a jet in the standard 
model is reflected by a statement in the SM constructor: 

construct _f ami ly ( " j " , " u , ubar , d , dbar , c , cbar , s , sbar , b , bbar , g " ) ; 

making the j name available in multi-process algorithm instances like CM_algorithin<SM,2,2> 
("e+,e- > j , j"). Again, the loadO method will ensure only unique processes are inserted 
in the list. Secondly, a model inherits the decouple.particle () and decouple_vertex() 
functions. Because erasing a particle or vertex from the global instance's lists is rather unsafe 
(previously built algorithm instances become invalid), the decoupling functions only set a flag 
to the particle or vertex, omitting their contributions from future amplitude evaluations. 

3.5 Implemented models 

There are several predefined models in the package. Some toy models for checking purposes 
are phi3, phi34 and WZM, constituting respectively a c/)^- and {(j)^ + i;^'')-theory and the Wess- 
Zumino model [8J. More physically relevant models are 

• scalar _QED: scalar quantum electrodynamics with a single flavor of scalar electrons. 

• QED: contains three flavors of leptons, of which only the r is massive by default. 

• QCD: contains six quark flavors, of which the c, b and t are massive by default. The static 
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member function set_f our_gluon_vertex() allows to switch between an fom'-gluon 
vertex description of the Lagrangian and an auxiliary tensor field treatment. 



• EWSM: the colorless electroweak standard model implementation, with massive c, b, t and 
r and both ud- and 6c-mixing. 

• SM: the full standard model, combining all features of QCD and the EWSM classes. 

• susy_QED,susy_QCD: supersymmetric quantum electro- and chromodynamics. The spar- 
ticles have by default the same mass as their standard model partners and zero decay 
width. 

All of the above model classes contain the following public static members: 

• alpha, alpha.s and QCD_scale: electromagnetic and strong couplings, and the scale 
at which as is evaluated. The former is constant and equals -1 in (susy_)QCD, the 
latter two in the colorless models. The functions set_alpha(. . .) , set_alpha_s (. . .) , 
and set_QCD_scale (. . .) insert new values for the respective parameters and update 
vertex couplings depending on them. 

• ref resh.couplings () : calculates all couplings from the current values of the basic 
parameters. The basic parameters are the strong and electromagnetic couplings, and 
for the EWSM and SM classes the additional input values G_F (Fermi's constant), M_Z (Z 
pole mass) and M_hO (Higgs pole mass) using the conventions of [9j. 

• ref resh_masses : although the values of the masses are passed by reference and 
recursive relations therefore always use the current value of a parameter, switching from 
a zero mass to a nonzero value is correctly updated by the algorithm once this function 
has been called. This is because the choice between the massless or massive polarization 
routines is not checked at each amplitude evaluation. 

• set_unitary_gauge : replaces all massive and massless vector propagators respectively 
by 

p^ — ' p^ 

In the electroweak sector, the would-be Goldtone bosons are decoupled. 

• set_R_xi_gauge : replaces all vector and would-be Goldstone propagators respectively 
by 

-ig^"" + z(l - Op^p^Hp^ - CM2) i 

The (electroweak) standard model classes also provide built-in functions for computing (lead- 
ing order) decay widths of the Higgs, top and W and Z bosons. These formulas assume all 
fermions massless except for the top and bottom quarks. Inclusion of finite widths is imple- 
mented by the complex mass scheme |10j . absorbing the width in a complex mass through 

M = \Jm'^ - iMT , 
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and performing the computation of the electroweak parameters with Mz, M\y and instead 
of the real pole masses. Because such a substitution preserves all Ward identities, resonant 
diagrams will still constitute gauge-invariant contributions to the amplitude. 

To simplify extensions all of the above models have templated counterparts called QED.- 
'bz.se< derived model,value type>, S}'lJDase< derived model,value type>, etc. Deriving a new 
model class from one of these automatically inserts their particle and vertex content in the 
user-defined class. Contrary to the model base class, these carry 2 template arguments. The 
second is the numerical type of the couplings and masses, and should be equal to the value_- 
type defined in the derived model. 



4 Library details 
4.1 Polarizations 

For the construction of massless spinor wave functions, CAMORRA adopts the spinor techniques 
developed in |lli I12j . From the beam.direction value in the model class, the algorithm 
constructs two momenta ko and ki that fulfill 

^2 = , kf = -l, ko-ki = 

and for which it assumes no momenta will be parallel to. For the common choice 3 (beam di- 
rection along the 2:-axis), these will be chosen to be fco = (1, 1, 0, 0) and ki = (0, 0, 1, 0). Then, 
the Dirac_algebra_type is used to construct two statically defined basic spinors u±{ko), ful- 
filling 

u±{ko)u±{ko) = ^(1 =F 7^)1^0 > ^^+(^0) = ^iU-{ko) , 

The massless fermion wave functions can be directly constructed from the constant spinors 
above. To cope with Majorana particles however, the relative phase between spinor and 
antispinors has to be fixed by the charge conjugation matrix, u\{p) = Cv^{p) [13]. The reason 
for this feature of the spinor phase choice will be explained in the next section. Hence, the 
algorithm will construct the positive- helicity spinors with the basic spinors, and the negative- 
helicity ones with the explicit charge conjugation: 

u+ (p) = ■ pu- {ko), u_ (p) = Cu^ (p) . 

V2p • ko' 

For the construction of massless polarization vectors, the beam direction constant is used 
once more to construct a third massless momentum k2 not parallel to any momenta. For 
the standard choice beain_direction=3, the choice k2 = ki + (1,0,0,0) is made and massless 
helicity vectors are constructed from 

1 



^x'^iP) = ^/g^^ ^A(p)7^^A(fc2) 



The spacetime class also defines a function splitting a massive momentum p into two massless 
vectors pi and p2 whose spatial parts are (anti-)parallel to p. The massive vector helicity 
states are then evaluated from 

e/(p) = -^ux{pih^u,{p2) , e,^ = . 
V 2m rn 
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For the standard choices of Dirac algebra basis and beam direction the formulas above were 
evaluated symbolically using MAPLE and the resulting expressions coded by hand. For 
nonstandard bases of e.g. the Dirac algebra basis, the contractions have to be performed for 
each event and CAMORRA will produce helicity amplitudes at a lower rate. Finally there are the 
massive spinors, for which the construction policy is determined by the spin_vector_type of 
the model. In particular, this class implements a function constructing for a given massive 
momentum p the spacelike spin vector s{p) for which p ■ s = Q and = —1. The package 
contains three implementations, which carry the names introduced in |14j : 

KS.type : s^{p) 

helicity.type : s'^(p) 

polarised.type : s^{p) 

where i denotes the beam direction. The massive spinor wave functions are then constructed 
according to the standard formula 

(^ + m)(l±75^(p)) 
2V^o • [P + ms{p)) 

and the antiparticle spinors by v\{p) = Cu^{p). Again, the spinor formulas for the three 
choices above were explicitly coded, avoiding multiple contractions. 

4.2 Fermion flow reversal 

In previous paragraph, we emphasized that the particle and anti-particle spinor phases should 
be consistent with the definition of the charge conjugation matrix. This is because CAMDRRA 
computes scattering amplitudes containing Majorana currents by explicitly reversing fermion 
flows [15J wherever it encounters a discontinuity. The strategy for treating Majorana fermions 
consistently in the recursive algorithm [1^ can be summarized as follows: 

• External Majorana fermions are always treated as Dirac particles (never antiparticles) , 
i.e. with the helicity spinors u\{p) and u\{p) for resp. initial and final states. 

• At each vertex, the program assumes an input Majorana current to behave like a fermion, 
that is, a column spinor for outgoing and internal currents and a row spinor for incoming 
currents. For consistency CAMORRA therefore ensures that at any vertex a produced 
Majorana current behaves like a fermion as well. 

• Whenever a fermion-flow discontinuity is encountered, the program reverses the Ma- 
jorana current by explicit charge conjugation. Therefore, models with no Majorana 
particles will never induce such reversals. When two combined Majorana currents give 
a fermion-flow violation, the first one will be reversed. 



p^ 

m p • /cq ' 
m|p| ' 
mip^ + m) 
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Figure 2: Recursive relations for vertices containing one Dirac fermion (denoted ■0 or ijj) and one 
Majorana fermion (denoted A). The charged bosons are denoted by phi^ and the vertex Feynman 
rule F is an arbitrary combination of Dirac matrices. The matrix F' is corresponding charge-conjugate 
matrix [TS], which acts on the Clifford algebra basis by (7'')' = ^7^) (7^)' = il^l^Y = 7^7^ and 



Just like the Fermi minus signs, the spinor reversal flags are determined at vertex tree con- 
struction. Note that a spinor reversal should not modify the current itself, but a copy of the 
subamplitude, since other vertex combinations-contributing to different Feynman graphs- 
might not require the current to be reversed. Therefore it was feasible to absorb the multipli- 
cation with the charge conjugation into the vertex contraction, resulting in recursive relations 
in e.g. fig. [2| 

4.3 Color flow optimization 

As already mentioned in section 2, CAMORRA accepts both adjoint and color-flow gluons and can 
run in either continuous or discrete color mode. In the adjoint color treatment, the recursive 
relations are optimized by omitting contractions over zero- valued color components, result- 
ing in a performance comparable to the color_f low mode. For discrete colors however, the 
continuous color flows can be traced through the vertex tree, and contractions over color ten- 
sors may be replaced with color flow combination rules [17] , often referred to as color-dressed 
recursive relations [18]. This increases the processing speed by an expected factor around 9, 
compared to the continuous case (see fig. [3|. Note that the recursive relations contain all 
explicit terms decoupling the (unphysical) trace part of the color-flow gluon currents 

= '^apGa (1) 

such that the resulting amplitude is exact to all orders in l/N^. The recursive relations from 
the available color structures in terms of the basis above were derived using FORM |19j , and 
are included in the package. 
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Figure 3: Amplitude evaluation times (in s) for n gluon amplitudes. Performance test was carried out 
by a 2.4 GHz Intel core 2 duo processor. 

4.4 Built-in generators 

Included in the library are basic Monte Carlo generators for the color, helicity and momen- 
tum degrees of freedom. For the phase space generation, the standard RAMBO algorithm 
|20j for uniformly distributed massless or massive momenta are available; for example, given 
a CM_algorithm object myalgo, one calls 

myalgo . set Jiomentum_generator<uiiif ormjiassless_ps , std : : random> (&Ecm) ; 

to generate massless momenta with a total invariant mass stored in Ecm. The second template 
parameter, std: : random represents a random number generator class, in this case a wrapper 
class for the standard library's rand function which we have defined in the std namespace. 
The pointer to the phase space generator can be recovered by calling the get jnomentuni_gen- 
erator member function. Cuts can then be applied as follows: 

using phase_space_variable : : eta; 

myalgo . get Jiomentum_gener at or ()-> insert _upper_cut<eta> (3 , feetamax) ; 

The code above restricts the third external particle momentum (the first outgoing one) to 
the r] < r^max range. It will not restrict the momentum generation itself, but rather throw a 
zero amplitude (without actual tree evaluation) whenever myalgo encounters an event which 
does not pass the phase space cuts. Phase space variables are accessed by the value function: 

double myeta=myalgo.getJttomentum_generator()->value<eta>(4) ; 

The helicity degrees of freedom are inserted similarly via the set_helicity_generator<i?/pe, 
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random number generator>0 member function. For helicities only the class unif orm_helic- 
ities is available for the type parameter. In the discrete helicity mode, this generator throws 
a dice, while in the continuous case it generates uniformly an angle (p to construct the wave 
functions 

e'^'^Uj^{p) + e~*'^ti_(p) , for spin-1/2 particles, 
e^^e+{p) + e~*'^e_(p) , for massless spin-1 particles, 
e"^e+(p) + Sq^p) + e~^'^£-{p) , for massive spin-1 particles. 

Finally, the color sampling is directed by the member function set_colour_geiierator<ij/pe, 
random number generator>0 . Again, the unif orm.colours type throws dice in the discrete 
color case and generates uniformly normalized complex color phases in the continuous case. 
More physically relevant are the adjoint.QCD and colour _flow_QCD classes. In continuous 
color models, the first generates real, normalized octets for external particles in the adjoint rep- 
resentation of SU{N) and complex normalized triplets for SU{N) fundamental-representation 
objects, whereas the second generates arbitrary SU{N) matrices for the gluons using eq. [!} 
In the discrete color mode however, the colour _flow_QCD type will generate constrained color 
configurations, conserving the number of colors and anti-colors of each type [TTj. The gen- 
eration proceeds by generating a color connection, i.e. an ordering of the anti-colors in the 
process with respect to the colors. After generating the colors uniformly, the anti-colors are 
determined by their position in the connection. The above color generation works for QCD- 
like theories, which have a constant static number N_c defined in its model class header that 
represents the number of QCD colors. 

4.5 Tests 

The core routines and implemented models of CAMORRA were extensively tested on an event- 
by-event basis. Because of the flexibility and control over the computation allowed by the 
program, many cross-checking has been performed. These can be performed by the user as 
well, by typing make check after building the library. The testing facility includes 

• Fermion contraction routine checks by comparing results in the Dirac and Weyl bases 
for gamma matrices. 

• Helicity state wave functions checks for normalization, mutual orthogonality and transver- 
sality. Comparisons between spin-summed amplitudes for different spin vector construc- 
tion policies. 

• Tree-building algorithm checks by counting diagrams and symmetry properties of am- 
plitudes. 

• Model implementation checks (especially for the EWSM and SM classes) by comparing 
amplitudes in different gauges. 

• Recursive relation consistency checks (especially of the Majorana Feynman rules) by 
rotating the choice of the final current in the vertex tree. 

• Color structure checks by comparing corresponding adjoint and color-flow amplitudes. 
Color combinations in the discrete color flow mode are checked against the continuous 
case. 
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• Direct checks of helicity amplitudes by alternative computations, such as and the 
vanishing Parke- Taylor amplitudes in QCD 121\ IH [22]. 

• Cross check between supersymmetric and ordinary QCD amplitudes [23] 




5 Conclusion and outlook 

We have written a scattering amplitude computation library that combines a full flexibility 
in terms of model implementation with the efficient recursive algorithm. This makes CAMORRA 
suitable for (tree-level) matrix element corrections to new physics processes with complex 
final states. Because the setup of the program is completely modular, new interaction and 
particle types may be included without modifying existing code. In terms of speed, CAMORRA 
is competitive with similar programs, whereas for multi-leg amplitudes, it exceeds model- 
independent programs based on the diagrammatic approach. Furthermore CAMORRA does not 
make use of external libraries (other than the standard C-|— |- one) or symbolic manipulation 
programs and allows a user with basic knowledge of the programming language to implement 
new models and compute helicity amplitudes. 

In the future, the library may be equipped with a suitable phase space integrator to yield a 
full, model-independent matrix element event generator. Such a device should simulate the 
peaking structure of the integrand efficiently to achieve a good accuracy in reasonable com- 
putation times. Recently it has been proposed that a backwards traversal of the vertex tree 
allows a mapping of the amplitude peaks onto phase space O [25| , using a multichannel at 
each vertex node. This would be the obvious strategy towards a generic phase space integrator 
for CAMORRA. 

Another extension of the library-which may be immediately implemented-is the construction 
of an MSSM model class, possibly with an interface to SUSY Les Houches Accord param- 
eters files \25\ I26j. In principle all the ingredients (color and Lorentz structures, Majorana 
fermions) are available and tested in CAMORRA, reducing the implementation to a translation 
of the Feynman rules [27] to the correct add.vertex statements in the MSSM class. 

Finally, interfaces with the FeynRules f28] and LanHep [29] packages may be feasible to 
compute helicity amplitudes for given Lagrangian densities, and which would also facilitate 
cross-checking the outcomes with other packages. This may require some new definitions of 
Lorentz or color structures and possibly the inclusion of higher-spin particles. 
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A The phase space interface 



The momentum, helicity and color degrees of freedom of an external particle in a process that 
is loaded and built by the algorithm class are contained in objects of type 

CM_algorit]im<model_type , N_in , N_out> : : phase_space_type 

whose addresses are accessed by the get_phase_space() member function of the algorithm 
class. References to the particle momentum components can for example be obtained as fol- 
lows: 

myalgo . get_phase_space (0) ->momeiituin(0) =450 ; 

This sets the time-component of the first particle's 4-momentum to 450 GeV. Input of helici- 
ties depends on the coiitinuous_helicities flag in the model class. If it set false, the phase 
space type contains a helicity integer that may be set as follows: 

myalgo . get_phase_space (2) ->helicity () =-1 ; 

Valid choices for the helicity are for scalars, ±1 for fermions and massless vectors and 
{—1,0, 1} for massive vector bosons. In the continuous helicity case, complex floating-point 
numbers can be assigned as coefficients belonging to the helicity states, e.g. 

myalgo. get_phase_space(2)->helicity_phase(-l)=std: :polar(3. 14) ; 

The color insertion proceeds in essentially the same way. If the continuous.colours flag 
is set false, the color integers are inserted by e.g. 

myalgo . get_phase_space (0) ->colour ( 1 ) =2 ; 

If the first particle in the process is a gluon, this line of code can be used for instance to 
set its anticolor to ^ (assuming the algorithm runs in color fiow mode). Note that particles 
with a single color index should always have it assigned by colour (0) = . . . (quarks do not 
have an anticolor in CAMORRA, there is no need to set it to zero). In the continuous color 
mode, the phase space instances contain color vectors or matrices (or tensors) of complex 
floating-point numbers. Their entries can be accessed by the function colour_coeff which 
takes a variadic argument list of integers: 

myalgo .get_phase_space(0) ->colour_coef f (1 ,2)=complex<'ya/ue type>(l , 0) ; 

is the continuous equivalent of the previous color assignment. Note that if the length of 
the argument list does not match the color rank of the fleld, the above code will result in a 
run-time error. 

^Remember that the colors run from to A'^c — 1- 
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B Examples 



The first example checks the vanishing Parke- Taylor amplitudes for a single 12-gluon event in 
QCD. The program instantiates an algorithm object called algo. After loading the process 
and constructing the vertex tree, it prints the number of Feynman diagram^ Then it uses 
the flat phase space generator at a center-of-mass energy of 1 TeV and the QCD color- 
flow sampler to generate momentum and color configurations. At this point all the gluon 
polarization vectors are initialized to e-{pi), and the amplitude is seen to vanish at double 
precision. It will remain zero after changing the wave function of gluon 2 to £-{p2) + £+(^2)- 
After including a second positive-helicity term, the amplitude becomes nontrivial. 

#include <Camorra/QCD . h> 
#include <Camorra/CM_algorithm.h> 
#include <Camorra/unif orni_inassless_ps .h> 
#include <Camorra/QCD_cols .h> 
#include <Camorra/rcarry .h> 

using namespace Camorra; 

int mainO 
{ 

CM_algorithm<QCD,2,10>algo("g,g > g,g,g,g,g,g,g,g,g,g") ; 
algo . loadO ; 
algo . construct ; 

std: : cout«"Nr of Feynman diagrams: "«algo.count_diagrams()«std: :endl; 
double Ecm=1000; 

algo . set_momentum_generator<unif orm_massless_ps ,rcarry> (&Ecm) ; 
algo . set_colour_generator<colour_f low_QCD,rcarry> () ; 
algo . generate ( ) ; 

std: :cout«algo. evaluate ()<<std: :endl; 

algo .get _phase_space (2)->helicity_phase(l)=l . ; 
algo .get_phase_space (2)->helicity_phase(-l)=l .0; 
std: :cout«algo. evaluate ()<<std: :endl; 

algo .get _phase_space (3)->helicity_phase(l)=l . ; 
algo .get_phase_space (3)->helicity_phase(-l)=l .0; 
std: :cout«algo. evaluate ()<<std: :endl; 

} 

After compilation the program will produce the output 

Nr of Feynman diagrams: 50363463150 
(-1 . 2957e-24 , -4 . 11328e-24) 
(-9 . 30578e-25 , 1 . 16322e-24) 

^Note that the 4-gluon vertex in the QCD class is replaced by an auxiliary antisymmetric tensor field 
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(-2 . 01376e-15 , -6 . 89148e-14) 



Note that to ensure portability we used the rcarry random number generator [30]. For 
completeness, the generated colors and momenta producing this output are included in the 
package. 

The second example involves the implementation of new model class. The toy model is called 
SMseesaw and adds a TeV-scale Majorana electron-neutrino nu_e to the standard model 
class. The model header SMseesaw. h file reads 

#ifndef SMSEESAW_H_ 
#define SMSEESAW_H_ 

ttinclude <Camorra/SM_base .h> 
#include <Camorra/Minkowski .h> 
#include <Camorra/Pauli_basis .h> 
#include <Camorra/helicity_type .h> 
#include <Camorra/col_f low.h> 

using namespace Camorra; 

class SMseesaw: public SM_base<SMseesaw,double> 

public : 

typedef double value_type; 
typedef Minkowski_type spacetime_type ; 
typedef Pauli_basis Dirac_algebra_type ; 
typedef helicity_type spin_vector_type ; 
typedef colour_flow colour_treatment ; 

static const unsigned dimension=4; 

static const unsigned N_c=3; 

static const int beam_direction=3; 

static const bool coloured=true ; 

static const bool continuous_helicities=f alse ; 

static const bool continuous_colours=f alse ; 

static value_type M_N; 

SMseesawO ; 

}; 

#endif 

The new model class is derived from the SM_base class to load all the standard-model particles 
before adding new physics. In the header declaration we see the usual compile-time type 
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definitions and constant static enumerable values. Note the N_c integer, which is necessary to 
define when deriving from QCD_base, SM_base or susy_QCD_base as it determines the number 
of colors in the QCD sector. Note also the static M_N floating-point number, which will store 
Majorana neutrino mass in the source file SMseesaw. cpp, 

#include " seesaw. h" 

using namespace Camorra; 

const unsigned SMseesaw: : dimension; 

const unsigned SMseesaw: :N_c; 

const int SMseesaw: :beam_direction; 

const bool SMseesaw: : coloured; 

const bool SMseesaw: :continuous_helicities; 

const bool SMseesaw: : continuous_colours ; 

SMseesaw: : value _type SMseesaw: :M_N=1000; 

SMseesaw: :SMseesaw() 
{ 

add_f ermion("~nu_e" ,&M_N) ; 

const std: : complex<double>*coupling=&SM_base<SMseesaw,double> : :Wnee; 

add_vertex<vffR>("W+" , "~nu_e" , "e-" , coupling) ; 
add_vertex<vffR>("W-" , "e+" , "~nu_e" , coupling) ; 

} 

The file begins with the initializations of the static data members of the model class. In 
the constructor, it adds a color-singlet Majorana electron-neutrino with mass M_N and width 
zero. Then it constructs a vertex between the electron, W bosons and Majorana neutrino, 
with couplings equal to the VKz/^£-couplings in the base class. After compiling this class 
into an object file, one can compute its predicted scattering amplitudes, for example for the 
d,u u, d, e~,e~ process: 

#include " SMseesaw. h" 
#include <Camorra/CM_algorithm.h> 
#include <Camorra/unif orm_massless_ps .h> 
#include <Camorra/unif oriii_hels .h> 
#include <Camorra/QCD_cols .h> 
#include <Camorra/rcarry .h> 

using namespace Camorra; 

int mainO 
{ 

CM_algorithm<SMseesaw,2,4>algo("d,ubar > u,dbar,e-,e-") ; 
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algo . loadO ; 
algo . construct ; 

std: : cout«"Nr of Feynman diagrams: "«algo. count .diagrams ()«std: :endl; 
double Ecm=500; 

algo . set_momentum_generator<unif orm_massless_ps ,rcarry>(&Ecm) ; 
algo . set_colour_generator<colour_f low_QCD,rcarry> ; 
algo . generate () ; 
algo . sum_spins () ; 

std: :cout«algo.evaluate_spin_sum()«std: :endl; 

> 

which counts the number of diagrams and computes the spin-summed amphtude of random 
phase space point at a 500 GeV center-of-mass energy of the incoming quarks. It should yield 
the outcome 

Nr of Feynman diagrams : 4 
1.69203e-12 

The momentum and color configuration yielding this amplitude is included in the file example2 . 
to be inserted manually in case rcarry turns out to give different results on the user's sys- 
tem. The examples above and the phase space information files are included in the examples 
directory of the package. 

C Installation and system requirements 

CAMORRA can be built on UNIX-based systems that arc equipped with GNU make and libtool. 

It does not require any additional packages other than the standard library. The installation 
is straightforward: after downloading the archive to a directory of choice, unpack and decom- 
press it by typing 

tar -xzvf camorra-l.O.tar.gz 

in a terminal window. Then, move into CAMORRA's main directory and type 

./configure — pref ix=<msia/Z directory> 
make 

maJse check (optional, takes about Ih) 
maJse install 

cd examples; make (optional) 

where install directory denotes the folder where the headers and library will be copied to. 
The — prefix=. . . may be omitted by users with root privileges, in which case the library 
will be installed in /usr/local/lib and the headers in /usr/local/include. More detailed 
information kan be found in the README and INSTALL files. In the doc directory of the 
distribution one finds a html reference list of all relevant classes and functions, generated with 
doxygen. 
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